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Hyperfine anomalies in Fr: boundaries of the spherical single particle model 
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We have measured the hyperfine splitting of the 7P\/2 state at the 100 ppm level in Fr isotopes 
( 206 S ,206m,207,209,213,22i Fr) near the dosed neutron shel i (A r = 126 in 213 Fr). The measurements in five 
isotopes and a nuclear isomeric state of francium, combined with previous determinations of the 7Si/ 2 
splittings, reveal the spatial distribution of the nuclear magnetization, i.e. the Bohr-Weisskopf effect. We 
compare our results with a simple shell model consisting of unpaired single valence nucleons orbiting a 
spherical nucleus, and find good agreement over a range of neutron-deficient isotopes ( 20 '~ 213 Fr). Also, 
we find near-constant proton anomalies for several even-A isotopes. This identifies a set of Fr isotopes 
whose nuclear structure can be understood well enough for the extraction of weak interaction parameters 
from parity non-conservation studies. 

PACS numbers: 21.10.Gv,27.80.+w,32.10.Fn 


Weak interaction studies in heavy atoms require for their 
interpretation precise knowledge of the atomic and nuclear 
wavefunctions. To extract nucleon-nucleon weak interac¬ 
tion couplings from the weak interaction induced parity- 
violating anapole moment (T), nuclei with simple and reg¬ 
ular magnetic properties are desirable [2]-|4). The nuclear 
magnetic moment is used to benchmark nuclear structure 
theories for calculating the anapole moment [0, which is a 
contact field effect produced inside the finite extent of the 
nucleus. Here we explore the regularity of the magnetic 
properties of a chain of Fr isotopes and find that 207 213 Fr 
in the vicinity of the neutron shell closure mark a range 
where the nuclear structure is sufficiently tractable for stan¬ 
dard model tests and constraints on new physics 0. 

To lowest order, the atomic hyperfine interaction can 
be described using a point-like nucleus characterized by 
the magnetic dipole moment. Deviations from the point¬ 
like approximation of the nucleus, referred to as hyperfine 
anomalies, come from considering how finite magnetic and 
charge distributions affect the interaction between the mag¬ 
netization of the nucleus and the magnetic field created by 
the electrons. The magnetic contribution is known as the 
Bohr-Weisskopf (BW) effect |6|[7). The difference in the 
nuclear charge distribution (Breit-Rosenthal (BR) effect 
If23l - i25l ) produces very small variations between isotopes, 
leaving the BW effect dominant Il26l l27il. As a new gen¬ 
eration of proposed parity violation experiments in atoms 
(including Fr) and molecules starts lf8l fl4l . it is important 
to understand the limiting factors due to the nuclear struc¬ 
ture, e.g. the nuclear magnetization, for the interpretation 


of parity-violating anapole moments |001- 

Here we present a systematic study of the variations of 
the hyperfine splittings (HFS) in a chain of Fr isotopes. Our 
measurements include 213 Fr with the closed neutron shell 
(magic number N = 126), 2069 ’ m Fr, 207 Fr, and 209 Fr on 
the neutron-deficient side, and the neutron-rich 221 Fr. In 
the region of 207-21,3 Fr with up to 6 neutron holes, we find 
near-constant magnetic hyperfine anomalies for the odd -Z, 
even-A 7 isotopes |fl5l . The neutron rich odd-even isotope 
221 Fr shows a different behavior due to the deformation of 
the nucleus. The odd -Z, odd-iV isotopes have anomaly 
contributions from both the proton and the valence neutron. 

BW effect measurements usually require precise knowl¬ 
edge of both, hyperfine structure constants and magnetic 
moments. Measurements of the nuclear magnetic dipole 
moment in Fr are limited to 211 Fr jT6l and 210 Fr ||T7I . Em¬ 
pirical values for other isotopes are usually obtained by 
scaling with the isotopic ratios of the hyperfine constants 
based on these two experiments, both of which have un¬ 
certainties larger than 1 % and cannot be used to extract the 
hyperfine anomaly. A different approach to study the BW 
effect which circumvents the limited precision in the mag¬ 
netic moments comes from the suggestion by Persson lfl8l 
that was implemented in Fr |fI51 , T1 ifTfil and Eu |[20) . This 
method relies on looking at the ratio between hyperfine 
splittings of two levels where the electron wavefunctions 
are overlapping differently with the nuclear wavefunction, 
and consequently is sensitive to the differences in their hy¬ 
perfine anomalies. We can normalize the change in this 
ratio to a specific isotope ( 213 Fr) for purpose of directly 
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revealing the contribution from the neutron wavefunction. 
See Ref. Q for a review and Ref. ET1 for a recent compi¬ 
lation of all available hyperfine anomaly data. 

The magnetic hyperfine interaction W can be written as 

DEED 

^extended = Wpoint(l + 6 i), (1) 

where ei is a small quantity that depends on the particular 
isotope, and on the atomic state (l = S or P). The 7 Pi /2 
electron overlaps with the nuclear wavefunctions more uni¬ 
formly than the 7 Si/ 2 electron. The ratio R of the hyper¬ 
fine splittings for an isotope with mass number A is sensi¬ 
tive to the nuclear magnetization distribution fl5l 

W s ^ , 

Rkfs(A) = - tended « R 0 (l + e s (A) - e P (A)),(2) 

extended 

with R 0 the ratio of hyperfine structure constants for a 
point nucleus. Since both states have J = the extraction 
of precise magnetic hyperfine structure constants from the 
measurements is not hampered by the presence of higher 
order nuclear moments. The relative size of e P grows with 
nuclear charge number Z , and is about 1/3 of e P in Fr l22l . 

We measure the 7P\ /•> hyperfine splitting of Fr at the 
100 ppm level in a number of isotopes. We use these mea¬ 
surements in combination with the 7Si/ 2 hyperfine split¬ 
tings mm to determine F?hfs(^ 4) to study changes in 
the hyperfine anomaly. These measurements are carried 
out at the Francium Trapping Facility (FTF) at TRIUMF 
fl8l . We briefly summarize the operation of the FTF: A 500 
MeV proton beam irradiates a target that consists of ura¬ 
nium carbide foils to produce between 10' to 10 9 Fr + /sec 
of the selected isotope ll28l . We produce an ultracold sam¬ 
ple of neutral Fr atoms for Doppler-free spectroscopy by 
capturing typically a few 10 5 atoms in a magneto-optical 
trap (MOT). Two Ti:Sapphire lasers (trap and repumper) 
form the MOT on the D 2 line (718 nm) and leave the D \ 
line (817 nm) background free for the measurement (Fig. 
|Tja)). A computer-controlled Fabry-Perot cavity monitors 
and stabilizes |j29l the long-term frequency variation of all 
of the lasers to better than ±5 MHz. 

A third Ti:Sapphire laser excites the transition from the 
upper ground state hyperfine level to each of the two 7P\/ 2 
hyperfine levels at 817 nm (Fig. |T|a)). We amplitude- 
modulate the probe laser with a fiber electro-optic modu¬ 
lator (EOM, EOSpace AZ-2K1-10-FPA-FPA-800-UL) that 
suppresses the carrier and produces sidebands at about half 
the hyperfine splitting, such that the two sidebands are sep¬ 
arated by the size of the splitting l30l . With this method 
we produce rapid scans of the sidebands OTj that min¬ 
imize many systematic effects compared to scanning the 
carrier frequency lfl5l . For instance, we are less sensitive 
to atom number fluctuations or laser frequency drifts, and 
the frequency axis of the scan can be precisely character¬ 
ized since it lies in the microwave regime. We produce 
the microwaves with phase-locked-loop synthesizers refer¬ 
enced to a Rb clock (SRS FS275). The frequency sweep 
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Time bin LJ 640 ns 


FIG. 1. Measurement scheme, (a) Atomic energy levels relevant 
for trapping and measuring, (b) Sideband frequency scan around 
2.9 GHz (isotope-dependent), (c) Time sequence for trapping and 
data collection. 


covers 140 MHz in 10 ms (Fig. |TJb)). The probe beam 
has a 3 mm diameter with 100 /iW power in each sideband 
and is retro-reflected to minimize trap displacement from 
radiation pressure imbalance. 

We collect the fluorescence light with a double relay 
imaging system (numerical aperture of 0.12) with an in¬ 
terference filter at 817 nm and an edge filter at 795 nm to 
suppress background light, in particular the trapping laser 
at 718 nm. A photomultiplier tube (Hamamatsu H7422- 
50) operating in photon counting mode detects the fluo¬ 
rescence, and we record the photon events as a function 
of time with a multichannel scaler (SRS SR430). We use 
a typical bin width of 640 ns with count rates below 250 
kilocounts/s. We avoid the ac Stark shift from the MOT 
trapping light with an experimental cycle of 27 /i s of trap¬ 
ping followed by 2.97 //s of probing with the trap laser 
off (Fig. |1}c)). Data with signal-to-noise ratios of > 20 
are obtained within a few seconds. Fig. [2] shows a typical 
spectrum, which yields a HFS splitting with statistical un¬ 
certainty at the 30 kHz level. The two peaks indicate the 
modulation frequency where the +1 (-1) sideband is reso¬ 
nant with the upper (lower) hyperfine peak. 

The linear Zeeman effect is the dominant systematic 
uncertainty, because the quadrupole magnetic field of the 
MOT stays on during the measurement, and the atoms pop- 
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Frequency (GHz) 

FIG. 2. Sample of experimental data. Top: Fluorescence counts 
as a function of sideband frequency for 213 Fr with separate 
Lorentzian fits to each of the resonances. The hyperfine splitting 
is the difference of the frequency position of the two peaks and 
it gives 5739.60 ± 0.030 MHz. Bottom: Normalized residuals of 
the fits. 


ulate all Zeeman sub-levels. We characterized this effect by 
changing the magnetic field gradient (7 to 15 G/cm), probe 
laser polarization (linear as used in the measurement, to 
circular), and the position of the atom cloud (1 mm dislo¬ 
cation). We put an upper bound of 540 kHz systematic un¬ 
certainty from the Zeeman shift. This agrees well with aux¬ 
iliary tests in Rb. The sideband laser power changes during 
the scan due to RF power variations, which contributes to 
a background structure in the data. We take scans without 
the trapped atoms to record this background. These varia¬ 
tions, in addition to line-asymmetries caused by laser drifts 
during the data accumulation, create small line-shape de¬ 
viations from a Lorentzian function, as can be seen in the 
residual structures of the fit in Fig [2] We evaluate a system¬ 
atic uncertainty of 100 kHz from the line-shape distortions. 
The differential ac Stark shift from the trap laser is miti¬ 
gated by the fast chopping technique, and the shift from 
the repump laser is 90 kHz with an uncertainty of 60 kHz 
(32) . Other potential systematic effects include Doppler 
shifts, probe laser power, frequency calibration and linear¬ 
ity of the scan, which we all evaluated to be at a negligible 
level. The total systematical uncertainty is 552 kHz. 

Table [I] shows the hyperfine splittings of the 7Pi/2 state 
for the five isotopes from this work, together with those 
reported in Ref. fl5) . The uncertainty includes both the 
statistical and systematic error contributions stated above. 
For 2t)ll Fr and 221 Fr we find good agreement with Ref. 
1151 and Ref. ||33l . respectively, with smaller error bars. 
The table also lists the ratio 7 ? H fs(^ 4) introduced in Eq. 
[2] based on the literature values for the 7S\ /■> splittings 
1341437) . The normalized ratio of the hyperfine anomalies 


TABLE I. Isotope, spin (/), 7P\/ 2 hyperfine splittings (HFS) and 
ratio 7 ?hfs( j 4) °f 7<S' 1 / 2 and 7P 1 / 2 splittings for 206 - 213 ' 221 p n 
We illustrate the neutron orbital configuration (v orbital) and its 
spin alignment with respect to the total nuclear spin (v spin). 


Isotope 

I 

HFS(7P 1/2 ) 

Rhfs(^4) 

v orbital 

v spin 

206 9 

3 

6009.14(55) 

7.6022(14)" 

a P3/2 + & / 5/2 

4 v7i 

206 m 

7 

6521.56(57) 

7.6086(10) c 

fb/2 

t Mr 

207 

4.5 

5559.04(55) 

7.6308(12) d 



208 

7 

6561.0(2.3)“ 

7.6053(30) c,d 

fb/2 

Mi 

209 

4.5 

5639.5(1.0)“ 

7.6307(16) d 



209 

4.5 

5638.36(56) 

7.6323(13) d 



210 

6 

6150.9(1.3)“ 

7.6035(17) d 

h/2 

Mi 

211 

4.5 

5710.5(1.0)“ 

7.6297(15) d 



212 

5 

6556.0(1.0)“ 

7.6042(17) d 

Pl/2 

Mi 

213 

4.5 

5739.43(55) 

7.6292(18) e 



221 

2.5 

2431.0(55) 

7.6581(26) d 



221 

2.5 

2433.0(3.9)^ 

7.652(12) d 




“ Ref. (131 . The 7,S'| j 2 values come from b Ref. lf34l . 0 Ref. 
(351 . d Ref. £36), e Ref. (37). For 208 Fr, the HFS(7S 1/2 ) is the 


weighted average and scaled uncertainty of the two 
measurements following Ref. 1381 . 2 Ref. (33). 

/(nFs(-4 )//i > hfs( 2I3) (with 213 Fr taken as the reference 
isotope for convenience), in a chain of isotopes A = 206- 
213 and 221 is shown in Fig. [3] The isotopes span neu¬ 
tron numbers between N = 119 to 134. The red squares 
correspond to measurements from Ref. lfl~5) . and the blue 
diamonds are the new results. For A = 206 we measured 
both the low spin (/ = 3) nuclear ground state 206,/ Fr, and 
the first high-spin (/ = 7), long-lived isomeric state 206m Fr 
(lifetime>10 s, deduced from MOT lifetime). We clearly 
distinguish 206m Fr from 2(l6fl Fr in the MOT due to their dif¬ 
ferent trapping and repumping laser frequencies. 

The shell model explains reasonably well the magnetic 
moments of the light Fr isotopes near N = 126 09) . We 
consider only the dominant orbitals of single nucleons in 
the shell model to calculate hyperfine anomaly differences 
for the 75 i / 2 and 7P 3 / 2 electronic states 1171 122) . The total 
anomaly combining the proton and the neutron is given by 

m 

e S ~ e P = e 7r/?7r + t-vPv, (3) 

where /3 Vt „ are the fractional contributions to the magnetic 
moment from the proton and the neutron, respectively. The 
calculated anomalies (in %) for a valence proton (e n ) or 
neutron (e„) for nh 9/2 , vpi/ 2 , vf 5/ 2 , and vp 3/2 orbitals 
are —0.57, —3.13, —2.75, and —1.75 respectively. The 
total anomaly is represented by the green circles in Fig. [3] 
Even though the single particle neutron anomalies are a few 
times larger than that of the proton, the neutron fractional 
contribution 0 V is typically only 15%, depending on the 
orbitals, yielding a contribution from the valence neutron 
of about the same size as the calculated proton anomaly. 
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FIG. 3. (color online) _Rhfs(^) as defined in Eq. 2, for sev¬ 
eral Fr isotopes normalized to the value for 213 Fr (closed neutron 
shell N=126). Blue diamonds are current measurements. Red 
squares are from Ref. G5J. The green circles are predictions 
based on Ref. lf22l . The ground state 206s Fr is shown to the left 
of the isomer 206m Fr. Both are enclosed in rectangles with the re¬ 
spective calculations. The proton resides on the indicated irh 9 / 2 
orbital. 


Our data shows that the even-odd staggering trend 
is preserved from 213 Fr with a closed neutron shell, 
down to 2ll< ’Fr. The neutron-deficient odd-even isotopes 
207,209,211,2 isp r a || have near-constant anomalies, differ¬ 
ing by a small amount consistent with variations due to 
the changes in the charge distribution (BR effect) [231127I . 
This confirms the previous assumption that the anomaly in 
these isotopes is due primarily to the single valence pro¬ 
ton in the nhg /2 orbital of the spherical shell [fl5l . 221 Fr is 
also odd-even, with the valence proton in the same irh 9 / 2 
orbital. Flowever the 8 neutrons above closed shell create 
deformation, leading to a rather different anomaly, which 
we explain in more detail later. For odd-odd isotopes 
( 208 , 2 io, 2 i 2 p r ^ t p e ca i cu i a tjons for the respective neutron 
orbitals are in good agreement with the experimental data. 

For 20 flfl Fr, there are some differences compared to a 
pure zzp 3/2 orbital, which also holds for the magnetic mo¬ 
ment fl34] i35 ., 40’l. Flowever other nearby orbitals (zz / 5 / 2 or 
7 /ii 3 / 2 ) give very different values, or even opposite signs 
with respect to the normalization. The calculations using 


different orbitals for 206 ®Fr are shown in Fig. [|] We note 
that the odd-odd isotopes have the same sign and roughly 
the same value of the neutron anomaly contribution. This is 
a coincidence resulting from the angular momentum cou¬ 
pling, even though their neutron orbitals and nuclear spins 
are different. This is illustrated by the fact the calcula¬ 
tion for lower nuclear spin (/ = 3) 12 / 5/2 ( 206 s Fr) pro¬ 
duces an opposite sign compared to the higher spin iso¬ 
topes ( 2 06m,208, 2 iop r ^ as / = 3 demands the neutron spin 
to be anti-aligned with the total nuclear spin (see the last 
column of Table [I). 

Deviating from the spherical shell model, the Nilsson 
picture fill considers the nuclear energy level changes due 
to a deformation. Calculations [43ll have shown that the 
f 5 / 2 and p 3 / 2 neutron orbitals have a level-crossing of the 
[521 1/2] and [503 5/2] sub-states at a very small (nega¬ 
tive) deformation parameter |e 2 | < 0.05, a range consis¬ 
tent with the measured quadrupole moment 051 . The va¬ 
lence neutron of 206 Fr can have a mixture of these two or¬ 
bitals, which shows effects in the magnetic moment [35). 
The anomaly is more sensitive to this orbital mixture, and 
could have a value outside of the range delimited by the 
two pure orbital predictions. This is because operator eval¬ 
uations such as (/5/2 |t" 2 1.P3/2) break the orthogonality of 
the two eigenfunctions, resulting a linear dependence of 
the anomaly to such mixtures, compared to the quadratic 
dependence of the magnetic moment. 

The neutron-rich odd-even isotope 221 Fr has only the 
proton anomaly contribution. Its Nilsson deformation pa¬ 
rameter lies in a positive small to intermediate range, 0.1 < 
e 2 < 0.2 fll6l . and the valence proton still occupies the 
h 9 /2 orbital. However, its angular momentum projects on 
to the nuclear symmetry axis, such that the nuclear spin 
becomes I = 5/2 [43]. Calculations in [43l yield more 
than 95% of the wavefunction in the zr/ig / 2 [523 5/2] state. 
Treating 221 Fr with this predominant contribution we ob¬ 
tain a correct sign of the anomaly with respect to 21 ? Fr, as 
shown in Fig. [3] and a magnetic moment of +1.85 // y, in 
reasonable agreement with the value +1.58(3) //.y empir¬ 
ically scaled from the hyperfine splitting lfl6l . We note 
that this is a considerably simpler picture than the nu¬ 
clear model used in Ref. HE). The Nilsson calculations 
are parameter-dependent, and the states with lower angu¬ 
lar momentum projections used in Ref. lTT6l would involve 
more contributions from configuration mixing effects. 

In conclusion we present precise measurements of the 
7 + 1/2 hyperfine splitting in several francium isotopes. The 
results allow us to study the hyperfine anomaly starting 
from a closed neutron shell ( 213 Fr N = 126) with a simple 
nuclear distribution, to the boundaries of the single-particle 
spherical shell model. We demonstrate high-quality spec¬ 
troscopic measurements, both with ground state nuclei as 
well as with an isomer. The present results provide the ba¬ 
sis for testing the validity and accuracy of future nuclear 
structure calculations, which will be necessary to extract 
weak interaction physics from parity non-conservation 
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measurements in francium. 
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